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Abstract 

We calculate the leading preasymptotic correction to the inclusive width 
b — > ccs (two massive quarks in the final state) due to the operator aH. It 
is found that this correction tends to cancel the 1/N C part of the inclusive 
width calculated using naive factorization. The absolute value of the effect 
is of order 0.25. We also discuss how quark-hadron duality is realized in the 
inclusive decays. 
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1 Introduction 



It was discovered recently |], that the dominant preasymptotic correction 
in the inclusive weak nonleptonic decays of heavy flavors is due to the opera- 
tor of the type aH || . It is of order mg 2 relative to the asymptotic (parton) 
expression (itiq is the mass of the decaying heavy quark) and occurs in the 
interference part of the amplitude. The sign of this correction is such that 
generically it tends to suppress the parton result for the interference part. 

In this note we further investigate the effects of the soft-gluon emission 
reducible to the operator oH . Two issues are considered: the transitions 
of the type b — > ccs, with two massive quarks in the final state, and the 
distribution in the invariant mass of the final quark- ant iquark pair. The 
analysis of the latter distribution allows us to trace, at least at the qualitative 
level, the origin of the leading preasymptotic correction in the interference 
part in terms of the exclusive channels |2|, [§, and relate it to the rule of 
discarding 1/N C in the weak exclusive decays ||, ^ . 

All notations and conventions are the same as in the first part of this work 
[JJ], hereafter referred to as I. The reader can also find in I an introduction to 
the preasymptotic corrections in weak inclusive decays and a representative 
list of references. 



2 The leading preasymptotic correction in b — > 

CCS 

As it was explained in I, the effective weak Lagrangian we work with has the 
generic form 

L eff = ^V^MBT^^a) + c 2 (5 i r^)(&r'W)}. (1) 

The constant Gf in eq. ([[]) is the Fermi coupling constant, Vi^ are the 
relevant CKM matrix elements, will be treated as numerical coefficients, 
= 7 M (1 + 75). In the transition at hand we identify 

A = b, B = c, p = s, a = c. (2) 

It is convenient to write the operators 0\ and O2 in the following form: 

o 2 = Ggr„A)(Bi» = (si»(ci», (3) 
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01 = ^0* + (4) 



1 = 2{(5V lx t a A){BV lx t a a) = 2(sT l ,t a b)(cT IM t a c), (Trt a t b = (l/2)5 ab ). (5) 

In the parton model the interference part of T(b — > ccs) is well known in the 
literature (see e.g. ref. ||). 



rffi = — - — G 2 F | V 1 V 2 | 2 Cl c 2 Im < H A j | / d 4 xiT{O 2 {x)O+{0)} \ H A ^ > 

(6) 



G2p 1 V X V 2 | 2 c lC2 M 5 F 3 (m7M 2 ) 



96vr 3 

where m is the mass of the c quarks while M is that of b. The subscript "int" 
marks the interference part of the amplitude. The function F% reflects the 
effect of the nonvanishing c-quark mass - it would be 1 if m=0. Actually, 

F 3 (x) = u (l - Ux - 2x 2 - 12x 3 ) + 24x 2 (l - x 2 )ln((l + v„)/(l - v )), (7) 

Vq = y/l - Ax, x = m 2 /M 2 . 

Observe that, in accordance with the general rule, r int ~ N® while the non- 
interference part is proportional to N^. 

The asymptotic (parton) expression for r int given in eq. (^J) obviously 
corresponds to the correlation function of 2 with the trivial component of 
Oi, namely that where the cc pair is in the color singlet state (see the de- 
composition in eq. (Q); in the parton model only the term 2 /N c is relevant). 
The leading preasymptotic correction is due to the second component of Of, 
namely 0\, with the cc pair in the color-octet state. This implies a soft 
gluon emission giving rise to the operator bT^G^t a b, see Fig.l. Calculation 
of the coefficient in front of the operator bT^G^t a b is rather straightforward, 
we comment here only on a few basic points. Some details relevant to the 
discussion of the invariant mass distribution are given in Appendix. 

Technically it is very convenient to carry out all computations using the 
background field method ( a pedagogical review is presented in ref. [0). We 
expand the quark Green functions in the background field keeping only the 
term 0(G°) and 0(G l ) where G is a generic notation for the gluon field 
strength tensor. Terms of higher order in G or terms with derivatives are 
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irrelevant to the operator at hand. The quark propagator in the coordinate 
space has the form 



O^TTi ,J*C-^ (Til * ^ 

We use the Schwinger gauge condition on the background field: 



{ / } + (terms with the even number of 7 matrices). 



x,A^x) = 0, Ap = gAlf. (9) 

The expressions in braces account for the quark mass; they reduce to 1 in 
the limit m — > 0; K\$ stand for the corresponding McDonald functions. The 
terms in S(x, 0) containing even number of 7 matrices drop out because of 
the chiral structure of the operators 0\ 2 ( T's in the vertices ensure that the 
only relevant part of S(x, 0) is the one with the odd number of 7 matrices). 

Now, using the general structure of the Green function it is very easy 
to show that the gluon emission occurs only from the antiquark line (Fig.l). 
The emission from the c and s quark lines vanishes. This property could have 
been foreseen on general grounds. Indeed, the operators 0\ i2 written in the 
form (§), force the gluon emission inside the cc loop. On the other hand, one 
could Fierz rearrange these operators prior to explicit computation. Then 
one will conclude that the gluon emission must occur inside the cs loop. 

The actual calculation is carried out in two stages. The correlation func- 
tion 

rffi = — - — G% I V 1 V 2 | 2 Cl c 2 Im < H A j I / d 4 xiT{O 2 (x)Ot(0)} \ H A ^ > 



(10) 

contains as a natural block the loop graph of Fig. 2 where the thick line 
denotes the c-quark Green function in the background gluon field, one of 
the vertices is V ^t" 1 while the other one is T v . Since we are interested only in 
terms linear in the position of the color matrix t a is unimportant; one can 
move it freely along the thick line. Let us denote by q is the total momentum 
of the cc pair. Once the imaginary part of the loop in Fig. 2 is found as a 
function of q 2 it is quite trivial to determine the complete imaginary part in 
eq. ([[(]) simply by integrating over the phase space of the s quark. 
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After these explanatory remarks let proceed directly to the results. The 
imaginary part of the Feynman graph on Fig. 2 is 

= + u)} (11) 

where 

s = q 2 , v = (l-Am 2 /s) 1/2 , (12) 

and D a is defined on Fig. 2. Notice that in the limit m — > the coefficient 
function in front of the braces in eq. (|TTD is actually 5(s); we will return to 
discussion of this point later. 

Furthermore, the imaginary part in eq. ([U]) reduces to that of D", 

ImzT{ / d 4 xO 2 (x)Ol(0)} = 2 J (Im^)^-^(^r^r^) (13) 

where p = {p°,p) is the on-mass-shell momentum of the s quark and the 
factor of 2 is the same factor that appears in the definition of the operator 
01 (see eq. ©). 

Integration over the phase space is readily performed, and we arrive at 
the following expression for T^: 

,(D _ C% | V X V 2 | 2 c lC2 M 3 . -~ . 5 



r = *' 1 Z ' 1 \. 3 < H A , | AG 0ll l l 5 A | H A , > F A (m 2 /M 2 



(14) 



where 



FAx) = v (l + (l/2)x + 3x 2 ) - 3x(l - 2x 2 ) In ^ + ^°^ (15) 

Ax 

and the notations are the same as in eq. (^). In passing to eq. (|14D we 
have restored the generic notations for the quarks and mesons, as in I. For 
instance denotes the meson with the quark content Ay, etc. It is worth 
reminding that A is the same as the b quark in the problem at hand. 

The function F A reflects the masses of c and c in the final state. If m=0 
F 4 (x = 0) = 1 and we reproduce the expression for Ar hadr from refs. [1], §]. 
As explained in ref. |3| the hadronic matrix element in eq. flUp (it is written 
here in the A rest frame, assuming M— > oo) is expressible in terms of the 
observable quantities 

< | AG 0ll l l 5 A | H M >= 2m 2 H M. (16) 
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Here m 2 aH = (3/4) (m^* — m 2 B ). Combining eqs. (||) and (|T6| ) we obtain 

r (o) M 2 F 3 (x)' 1 J 

At this point we pause to mention a subtlety associated with the parameter 
M — the quark mass. Unlike m 2 uH this parameter is not directly observable 
and, hence, requires comments. A constant shift in M, i.e. independent of 
the heavy quark mass, gives rise to a power correction analogous to the one 
in eq. (|T7|) , generally speaking. In other words one can say that redefining 
M 5 in eq. (|6]) one generates a universal power correction in the non-leptonic 
widths of all A-flavored hadrons. 

Since we are interested only in the leading power correction, the definition 
of the quark mass will not change r^. This is not the case in rj^, however. 
As a matter of fact, it is not difficult to see that the bare graph corresponding 
to r[°t (analog of Fig. 1, but with no soft gluon emission from the cc loop) 
is proportional to P A P where P is the momentum of the initial quark A. 
In terms of the covariant quantities one can identically rewrite |l[ (see also 
Appendix) 

AP A PA = A{V 5 - l -({G a pa a/3 )V 3 + P\G a ^ aP ))}A (18) 

where V is the covariant momentum operator, = iD^. It is instructive to 
notice that in the limit M —>■ oo the operator P a AG a pT 13 A appearing in eq. 
flU]) and the operator AG a pa al3 A in eq. (|I8"D reduce to each other. Namely, 
in the nonrelativistic limit 

AGoa^A -> -h A aHh Al AG aP a a(3 A -> 2ih A aHh A . (19) 



Eq. flT8|) implies (through equations of motion) that M 5 in eq. (§) is 
actually to be understood as Mq — 2MQm 2 rH where M is the "bare" quark 
mass entering equations of motion. 

In the above arguments the Schwinger gauge condition on the background 
gluon field is always implied. Had we chosen another gauge for the back- 
ground field additional graphs with the soft gluon line attached would appear. 
The final assertions are, of course, gauge independent. 

Concluding the section let us present a few numerical estimates. We 
assume the b quark mass M to be equal to 4.5 GeV and the charmed quark 
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mass m to 1.35 GeV. Then x ~ 0.1, F3 = 0.15, F4 = 0.23. For the ratio of 
the factorizable and nonfactorizable parts of the 1/N C piece of the amplitude 
we obtain 

r = = -8^f-^) J - ~ -0.25. 20 

m| F 3 (0.1) 1 ; 

It is seen that this ratio is numerically the same as the corresponding ratio 
in the b — > cud channel, as it was calculated in refs. [1], 0. 

3 Duality. The invariant mass distributions 

This section is devoted to related questions: (i) how the interference part 
of the inclusive width might be saturated by specific hadronic states; (ii) 
duality in the asymptotic regime and preasymptotic corrections. 

We will check the self-consistency of our approach. First, we will see that 
quark-hadron duality is valid at each step of our calculations. Second, we will 
demonstrate that the sum of exclusive channels nicely matches the inclusive 
width. 

The discussion of how duality might be implemented in r int was started in 
[H] (see also I). Here we present additional arguments which seemingly solve 
the problem at the qualitative level. The invariant mass distribution, to be 
considered below, has a remarkable structure, with far-going consequences. 

Let us recall that in our calculation of the nonfactorizable contributions 
to the exclusive widths we used the expansion in the parameters (see ||) 

x q*~> x= qi (21) 

where A is the energy release, /x is a typical off-shellness of quarks (gluons) 
inside hadrons and Q 2 ~ 1 GeV 2 can be considered roughly speaking as the 
invariant mass of the light quark pair forming the hadronic state at hand. 
If A' <C 1 and the operator product expansion is applicable we have two 
possible kinematical regimes for A: 

1) A<1, 



2) A > 1. 
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The first regime corresponds to exclusive decays analysed in [[3] while 
the second one corresponds to the inclusive description. Generally speaking, 
situation with the power corrections in these two domains is different, and 
one should always keep this difference in mind. Still, one might hope to get 
a sort of matching, and we analyse this possibility. In other words, it will 
be explicitly checked how our predictions for exclusive and inclusive widths 
match each other. 

Since the issue is quite general, to begin with it is convenient to get rid 
of inessential details and choose kinematics in such a way as to make all 
assertions most transparent. 

Thus, we consider a generic process A — > Bf3a. We consider the kine- 
matical regime A > 1, relevant to inclusive processes and try to understand 
how the inclusive rate is saturated by exclusive channels. For simplicity we 
impose an additional constraint (which helps technically but is otherwise 
unimportant), 

M A , M B -> oo, M A - M B = A fixed, m a =m p = 0. (22) 

This means that we to work in the small velocity (SV) limit |5j, avoiding 
complications due to recoil dependence of various formfactors. 

The central idea is as follows. In the asymptotic limit where the preasymp- 
totic corrections are negligible there exist a natural line of reasoning showing 
that neither the naive factorization nor the rule of discarding 1/N C can take 
place universally (see I). The actual situation is more complicated. Pre- 
sumably, for the operator 0\ = (BT fJ ,A)(PT ll a) the rule of discarding l/N c 
holds maximally while for the operator 0<j~ = (AT ll f3){aT^B) it is maximally 
violated, and the naive factorization takes place. [] The arguments substanti- 
ating the point are admittedly heuristic and are based on the notion of color 
transparency (see e.g. fllPp. 



Let us elucidate the above assertion in more detail. We consider now the 
limit of infinitely large Ma,b with the preasymptotic corrections switched 
off. (The parton- model limit). As explained in I the absorptive part of the 



1 We hasten to add that this asymptotic "inclusive" picture is quite opposite to what 
happens in the exclusive decay B° — > D + n~ . It has been shown || that the rule of dis- 
carding 1/N C is applicable to O2 in this decay. The breakdown of the expected asymptotic 
picture in the exclusive channels is correlated with the preasymptotic corrections to the 
inclusive rates, as will be discussed below. 
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interference amplitude is contributed, in the general case, by the intermediate 
hadronic states of 3 distinct types. 

(i) The operator d = (BV^A) (pVa) produces the state {£7} + {pa} 
in the leading order in 1/N C ; this state is absorbed by the operator = 
(AT ^P)(aT fl B) at the level 1/N C (the figure bracket denotes a color-singlet 
meson state, possibly excited, with the given quark content); 

(ii) the operator 0\ produces the state {Ba} + {P^} at the level l/N c 
which is then swallowed by O2 in the leading order. 

(iii) the operator 0\ produces the state {Pa} + {Bp} + {P^} at the level 
1/y/Nc] this state is then annihilated by 0% at the level l/\fW c . 

Let us argue that only the first option is realized in the kinematic domain 
considered. The state {Ba} + {^7} that formally has a projection on 0\ 
of order 1/N C actually has an intrinsic suppression by powers of Aq C d/A, 
where A is the energy release. (A = Ma — Mb). Indeed, in the SV point 
(vb — > 0) the quark B "gently" substitutes A. The {7-B} system remains 
undistorted. The a/3 pair produced at the time t=0 does not interact with 
the surrounding medium till t ~ Aqq D i n its res ^ frame, i.e. till t ~ A/Aq CD 
in the rest frame of A. By that time the a (3 system is separated from S7 
by distance of order A/Aq CD 3> A^ CD . From the point of view of the {/3a} 
system at this moment of time the state {-67} presents a tiny white object far 
beyond the correlation length. Therefore, it disappears "below" the horizon 
and decouples (Fig. 3). The color twisted state {Ba} + {f3^} does not evolve 
from 0\. The corresponding suppression is governed not only by l/N c ; there 
is another parameter, Aq C d/A. 

One can identically rewrite 1 using the Fierz transformation, 

Ox = ^T^A)(BT»a) + 2(PT,t a A)(Bm a a). (23) 

This means that at short times the operator may 1 form a superposition 
of two states of the type "/37" s + " Ba" s plus "/37" + "Sa" (where the 
subscripts s and o mark color singlet and octet), but at large times the octet 
x octet part cancels singlet x singlet , and the rule of discarding 1/N C holds 
for the production of color- twisted states by 0\ (Fig. 4). In other words, the 
states of the type (ii) do not play a role in the saturation of . 

Proceed now to the states of the type (iii). Here the situation seems to be 
the simplest. The production of the extra (3(3 pair from {-87} is suppressed 
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by l/y/Nc and | Vb \, where | vb | is the B-quark velocity in the A rest frame. 
Hence in the SV limit this contribution is absent. 

The last point to be discussed is the states of the type (i) and their 
projection on O^. The operator 2 after the Fierz rearrangement takes the 
form 

2 = ±-(BT^A)C^a) + 2(BT fl t a A)(PTH a a) (24) 

It produces the desired state at the level 1/N C plus the octet x octet com- 
ponent (Fig. 5). The latter includes a pair of fast light quarks. This latter 
component asymptotically vanishes, presumably due to the strong radiation 
of gluons (the Sudakov effect). Thus, for the color twisted states produced 
by O2 we have exact naive factorization at the level of 1/N C . (Warning: we 
speak here of highly excited states typical of the inclusive production, see the 
footnote on page 8). As a result of this asymmetry in the hadronic saturation 
(the rule of discarding 1/N C for 0\ and the naive factorization for O2) the 
duality in the interference part of the amplitude is now perfectly satisfied. 
The origin of the asymmetry between 0\ and O2 is specific kinematics, of 
course: the SV limit for B and large energy release - fast (3a. Under different 
kinematic conditions the pattern of saturation can be different ||. 

The naive factorization for the operator O2 that we advocate here is of 
course not in contradiction with the results of ref. ||, where the rule of 
discarding 1/N C has been shown to emerge for the same operator, but in the 
different kinematic regime. In ref. ||] the light quark pair we dealt with had a 
limited invariant mass, a few units times Aqcd- In the inclusive decays in the 
parton-model regime discussed above the distribution over the (3a invariant 
mass is of order A, and this parameter is assumed to be parametrically much 
larger than Aqqd- Therefore, the rule of discarding 1/N C for small invariant 
masses and the naive factorization for large invariant masses for one and the 
same operator may well coexist. Shortly, when we proceed to the discussion of 
the preasymptotic corrections, we will find additional arguments supporting 
this assertion. 

We saw above how duality is realized in r^. Let us now consider the 
preasymptotic corrections. 

Now we come to the most essential question— the invariant mass distribu- 
tion for the light quark pair (flat), produced by the operator 0\ and absorbed 
by O2. The point is that in the case of the soft gluon emission (operator aH) 
this distribution is drastically different from the one that occurs in the parton 
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part of the amplitude. Indeed, let us analyse the expression for IhlD^, given 
by eq. flTT|) in the limit m — > 0. Recall that s is the invariant mass squared 
of the (3a pair. For fixed s and m — > the parameter v defined in eq. (|T2]) 



tends to 1 and formally ImD* vanishes (see eq. flTT|)). It is not difficult to 
understand, however, that the above assertion of vanishing takes place only 
at s 7^ 0, and as a matter of fact at m — > ImD^ is proportional to S(s). 
Specifically, 

ImD^ \ m ^ = -^-5(s)(G ^ a q a q^ + (a* <-> ^)). (25) 

(Note the similarity of this expression to the Dolgov-Zakharov-'t Hooft ap- 
proach to the triangle anomaly JTT], [T3] ) . Eq. has a natural interpretation 
in terms of hadronic states. It implies that the hadronic states responsible 
for the leading preasymptotic correction are characterized by a small invari- 
ant mass of the (3a pair. Then, these quarks have to materialize in the 
form of light mesons, like 7r and p. In other words, the Dalitz plot, which is 
relatively homogeneously populated in the parton model is distorted in the 
domain corresponding to small invariant masses of (3a pair. The distortion 
corresponds to an extra negative contribution relative to that emerging in 
the factorized approximation. The operator O2 which obeys the naive factor- 
ization for highly excited states (populating the interior of the Dalitz plot) 
switches to the rule of discarding 1/N C for low invariant masses of (3a (i.e. 
at the boundary of the Dalitz plot). 

Indeed, if we consider the interference part of the amplitude, in the fac- 
torized approximation in the low s domain, the contribution of Dtt, Dp, etc. 
is obviously proportional to j^A . At the same time the preasymptotic cor- 
rection due to BR is ~ m^A 3 . The ratio of these two effects is N c m 2 aH / f%, 
i.e. exactly the same as has been found in ref. for exclusive decays of 
the type B — > Dti, B — ► Dp. Thus, the correction discussed in I and in 
the present paper can be most probably associated with the distortion of the 
decay modes where (3a pair forms generic it's and p's - dynamical emergence 
of the rule of discarding 1/N C in the exclusive decays of this kind. 

It is instructive to trace not only the parametric dependence like N c m 2 aH / f% 
mentioned above - which comes out correct - but to work out numbers as 
well. The contribution of Dir and D*n in the factorization approximation is 
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(neglecting the pion mass): 

Here and below we measure Ar in the units of G 2 ?\ViV2\ 2 ciC2- The Isgur- 
Wise function £ = 1 in the limit of eq. (22|). Moreover, if we neglect m p 



compared to A the contribution of Dp and D*p is the same as above, with 
an obvious substitution f 2 — > f 2 . According to ref. the nonfactorizable 
part of the same contribution for Dn and D*n is 

Ar „,._ m 3 (1 _g ) 3_l_ m J Hx(1+ l ) . (27) 

Here the factor 4/3 takes into account different degrees of cancellation of the 
1/N C parts in Dn and D*7r; we also omit all factors like x or x' which are of 
order one but are not exactly fixed by our calculation (see ref. ). For the 
case of the decay with the production of p the answer is essentially the same, 
up to corrections due to the nonvanishing mass of the p meson. Neglecting 
this mass we just double the result. Thus, we get, summing the contributions 
of these four channels in the nonfactorizable part of r int : 

Ar „, = - m%(l - ^)^ml H . (28) 

This result is remarkably close to the nonfactorizable contribution in the 
inclusive width rj^ calculated in []!], As a matter of fact, eq. (26) in I 
gives for r£jj the same expression with the substitution (1/24tt 3 ) -> (l/12vr 3 ). 
Thus, if we stick literally to our formulae in the exclusive case we find that 
the nonfactorizable correction in the four exclusive channels is responsible 
for ~ 0.5 of the inclusive nonfactorizable correction. The remaining 1/2 may 
be due to, say, a\ contribution. It may well be, however, that various uncer- 
tainties of order one in our exclusive estimates actually add up to increase 
the nonfactorizable parts of some of these four exclusive amplitudes, so that 
they build up the inclusive correction (almost) entirely. 

Summarizing, one can say that the nonfactorizable correction due to the 
operator bG^^Yl^b in the inclusive width nicely matches the corresponding 
sum in the exclusive channels. Quark-hadron duality is valid at each step of 
our calculation. 
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4 Conclusions 



In this paper we found that the leading preasymptotic correction to the 
inclusive width in the channel b — > ccs cancels about 25 % of the 1/N C 
suppressed part of the inclusive width calculated in the naive parton model. 
As in I, this effect leads to an additional increase of the inclusive hadronic 
width in this channel by ~ 5%. It would be interesting to have a direct 
experimental measurement of this width. 

We have also shown that the analysis of the invariant mass distribution 
allows one, at least at the qualitative level, to trace the way quark-hadron 
duality is realised in the weak decays of heavy quarks, and to understand 
the relation between inclusive and exclusive widths. In particular, we have 
demonstrated, as it could be expected already from the results of I that 
the leading preasymptotic correction to the inclusive width of B meson is 
saturated by the exclusive widths of the type B — > D + light meson, B — > D* 
+ light meson, where the set of light mesons seems to be exhausted by tt, p 
and may be A\ and p' . 

It would be of interest to calculate next-to-leading preasymptotic cor- 
rections both for b — > cud and b — > ccs channels and to study further the 
pertubative QCD corrections ]T3] . 

The authors are grateful to A. Dobrovolskaya, A. Vainshtein and M. Whit- 
terell for useful discussions. 



5 Appendix. Calculations in External Fields 
in the Interference Part of the Amplitude 

In this appendix we provide a few useful computational details relevant to 
the calculation of InxD" rj^ , and the covariantization of the parton model 

calculation of r^. 

We consider first the c-quark loop in the external field. This loop can be 
formally written as 

D a = zTr {iy a ^ T„-, } (29) 

1 ^ V-q/2-m V + q/2-m 1 1 ' 

Here V is the formal covariant momentum operator (see ref. 0). Since we 
are interested only in terms linear in G^, the gluon field strength tensor, 
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we can assume that G a v is constant (x-independent) c-number field. Using 
the general properties of the trace operation and the fact that in our ap- 
proximation t a can be moved to any position freely, it is trivial to show that 



Now, we can use the fact that 



(f> _ m )"l = (p + m )(p2 _ m 2 + l aapGaf} )-\ (30) 

The term with <J a pG a p is the magnetic emission, and this source of is 
quite obvious. One should notice, though, that G® is also hidden in V a \ this 
is the electric emission. Using eq. ( p0[) we rewrite eq. (|29|) in the following 
form 



1 

{V-q/2) 2 ~m 2 + UaG) 



D% = iTriT^iV - q/2) ^ ^ 2 _ 2 ■ if ^ T„ X (31) 



CP + 9/2) (p + q/2 y _ m2+ l _ (jG )- 

The mass terms drop out because of the chiral structure of V matrices. 

For the magnetic emission we merely expand the denominators of the 
quark propagators in aG keeping only the linear term in the expansion. 
Since we are aimed at linear terms only, in all other places one can immedi- 
ately substitute P — > p, the normal momentum, to be treated gular 
integration variable. Thus, calculation of this part is trivial. After some 
straightforward algebra we get in this way 

Im(Z^) magn = -M{G%q a q u + (// <-> u) + t(G%q a q„ + (// <-> i/))} x (32) 

1 3-v 2 
3g 2 ^ 2v 2 

Here Dq is the bare loop of scalar colorless particles: 



. r d 4 p 1 1 

°~~ t J (27r) 4 (p + g/2) 2 -m 2 X (p-q/2) 2 -m 2 ' ( ' 

ImD = -^, v = {l-Am 2 /q 2 ) 1 ' 2 . 

l07T 
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Situation with the charge emission is somewhat more complicated. Now 
we omit (crG) terms in eq. (|3~i"D. Then taking the trace over the Lorents (7 
matrix) indices and isolating the structure of interest we arrive at 



(iVW = 8iTr {f (P- g /2) _ _ m2 (7> + g/2) + _ w2 

(34) 

+ (1/ «-> /i)} 

It is worth emphasizing that eq. ( p4|) is automatically p-even (there is no 
p-odd contribution from the charge emission). It produces the structure 
(-D^) charge = -^ch (G° a q a q u + G a va q a q^). In order to determine the coefficient 
F c h in front of this structure it is convenient to apply the following trick. Let 
us introduce an arbitrary vector n a with the properties (nq) = 0, n a G a p 7^ 0. 
Then it is easy to obtain 

F ch (nG a q) = 4iTi{t a (Vn)x 

11 11 



(V - q/2) 2 -m 2 (V + q/2) 2 - m 2 (V + q/2) 2 - m 2 (V - q/2) 2 - m 2 



}• 

(35) 

In eq. (|35[) we omitted some tadpole pieces that have no imaginary parts. 
Let us take into account now that for any 2 operators 

kk - kk ^ kk°^ - °^kk- (36) 

This fact implies that in the square brackets in eq. (|35D we get the commu- 

[{{V - q/2) 2 - m 2 ), ((V + q/2) 2 - m 2 )] = 4iV a G a ^. (37) 

Now G a p is explicitly singled out, and we can substitute V — >• p everywhere. 
As a result 



F ch (nG a q) - 8G^y {2 ^P P P\ {p + q/2 y _ m 2y {{p _ q/2? _ 

(38) 

The integral in eq. (|38D can be done directly. Alternatively, one can reduce 
its imaginary part to IulDo by differentiating Dq with respect to q 2 and m 2 . 
Both derivatives are expressible in terms of d/dv. Moreover, since lmD is 
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linear in v the second derivative d 2 / dv 2 lmDo vanishes. In this way we arrive 
at 



(ImZ^) chargc = -Bi±(^^l v JL + ±)ImD (G; a q a q u + (/,<- u)) (39) 



= ~ 2l K 2 1) img {Gg n g a q, + (/* ~ ")}• 

In the sum corresponding to both magnetic plus electric emissions the p-even 
part cancels, and we are left with the p-odd part only, 

ImD^ = 2±^ImA>{G>% + (/i «-> i/)} (40) 



^7r q 2 v 



This expression, obviously, is due to parity violating combination of T's 

(7^75 X7, + 7 m x 7^75) in eq. 

It is instructive to notice that at m 2 —>■ InxD" , tends to a delta function 



ImZ^ | m2 _^-> ^(</ 2 ) + (// <- «/)}. (41) 

Let us proceed now to calculation of the full imaginary part in eq. (|10|). 
This calculation reduces now to (see Fig. 1) 

ImA4 = Im J d*xiT{O 2 (x)dt(0)} = (42) 
( ' '' ; lmD a (bt a T^pT u b) 



2p°(27r) 3 ^ y 

where p is the momentum of the massless s quark (p is on mass shell, p 2 = 0). 
The factor 2 in the r.h.s. reflects the factor of 2 in the definition of 0\. 

Integration over p is best done in the rest frame of b. Integrating over 
angles one finds 

where s = q 2 is the invariant mass of the cc pair, M is the b-quark mass. 
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Doing the integral in eq. ( ^3|) we arrive to eqs. (p!4[), (p!5|). 

It is instructive to check that the gluon emission from the line lying 
outside the closed quark loop vanishes, irrespectively of the relation between 
the masses of the quarks that propagate along the different lines. To this 
end it is sufficient to use eq. (§) for the quark Green function. The quark 
loop is formed from two bare propagators, while in the line lying outside the 
loop we take the 0(G) term in eq. (|]). The terms with the even number 
of 7 matrices (omitted in eq. (§) ) are irrelevant, since they drop out due 
to the chiral structure of T's. Using eq. (H) it is easy to see that the closed 
quark loop yields i^lx^Xy — x 2 g lll ) times a function of x 2 . The line lying 
outside the closed quark loop produces x a G aip bT u 'j lp 'y 5 T^b times a function 
of q 2 . Convoluting these two structures we get identical zero. 

Our last explanatory remark refers to eq. (|7|). Calculating the bare loop 
of Fig. 7 we get a local operator of the type bd 4 db, the intermediate quark line 
shrinks into a point. The operator bd^db contains derivatives which should 
be rewritten in terms of covariant quantities, say, for the purpose of using 
the equations of motion. 

Let us assume for definiteness that the operator bcftdb is taken at the 
origin (this assumption is technical and final result does not depend on it). 
The four-potential can be chosen in the form 

A; = \x?G pa . (44) 

The fact that A(0) = is the consequence of the Schwinger gauge condition 
imposed on the external field. Eq. ([14]) assumes this condition, and, if not 
stated to the contrary, all specific calculations made above are done in this 
gauge. In particular, the assertion that there is no A a ^ leg associated with 
the line lying outside the closed quark loop is also specific to this gauge. In 
the Schwinger gauge the result for the interference term < O2O1 > reduces 
to < O2O+ >~ Cr" plus the correlator jj- < O2, 0£ > given by the only 
graph depicted on Fig. 7. (modulo terms of higher order in A^) . Using eq. 

( f4*4"D and the remarks above it is easy to obtain the expression for P 4 P. We 
start from 

b(0)V 4 Vb{0) = 5(0)(P + Af{P + i)6(0), P = id^ (45) 
and commute all the A's to the left of all P's. Once A is in the leftmost 
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position we use the fact that A(0) = 0. In the Schwinger gau 



[PpA a ] - -Gf3 a . 



In this way we arrive at 



Besides that 



p4 fi = p 4p _ l[(Ga)P 3 - P 3 (Ga)]. 



V * = <p 2 - l -(Ga) 
Combining equations ([46] ) and fl47|) one finds 



P 4 P = V 5 - ^[(G(t)P 3 + P 3 (Ga)] + 0(G 2 ) 



which is equivalent to eq. (^). 
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Figure Captions 



Fig.l Imaginary part of such graphs determines the inclusive widths in 
the asymptotic regime. 

Fig. 2 A subgraph that determines the invariant mass distribution. 

Fig. 3 Creation of the states of the type (i) in the color transparency 
picture. 

Fig. 4 Creation of the states by the operator 0\ at the level 
Fig. 5 Creation of the states by the operator 2 at the level jj-. 
Fig. 6 Calculation of the correlator < 2 ,02 >■ 
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